Masking effect of heat dissipation on the current-voltage characteristics of a 
mesoscopic superconducting square with leads 
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A theoretical analysis based on a numerical solution of the coupled time-dependent Ginzburg- 
Landau and heat dissipation equations shows a strong dependence of the critical currents on the 
applied magnetic field in a mesoscopic square with attached contacts. In agreement with experiment 
we found hysteresis which are caused by a strong heat dissipation in the sample at currents close 
to the depairing Ginzburg-Landau current and /or the dynamics of the superconducting condensate. 
The theoretically obtained nonmonotonous dependence of the switching current (from superconduct- 
ing to the resistive state) on the applied magnetic field, arising from the changes in the vorticity, 
agrees quantitatively with the experimental data. Our results show that heat dissipation leads to an 
increase of the hysteresis in the current- voltage characteristic and hence masks the actual dynamics 
of the superconducting condensate. 
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I. INTRODUCTION 

After the discovery of the step-like features in 
the current-voltage characteristics of superconducting 
whiskersi and the explanation of this effect through the 
nucleation of phase slip centers 2 *^ this phenomenon at- 
tracted a lot of attention from several experimental and 
theoretical groups (see Refs4£ for reviews). This effect 
is a consequence of the nontrivial dynamics of the su- 
perconducting condensate and the normal quasiparticles 
near the phase slip centev^^. During the phase slip pro- 
cess the gap goes periodically in time to zero in one point 
along the sample and it creates an excess of quasiparti- 
cles (electron-like and hole like) near this region. Due to 
the relatively large time of relaxation of noncquilibrium 
quasiparticles in the superconductor they may diffuse on 
a distance much larger than the coherence length (the size 
of suppression of the gap)&. It results in the existence of 
a nonzero electrical field (as a response to the gradient of 
chemical potential of noncquilibrium quasiparticles) and 
a finite time-dependent voltage drop near the phase slip 
center. Unfortunately, strong heat dissipation masks this 
effect at temperatures far from the critical temperature 2 
T c which prohibited the study of this effect in full details. 

Recently, this subject was revisited because new exper- 
imental techniques were developed which made it pos- 
sible to prepare samples with low resistance^ (i.e. di- 
minishing the heating effects at low temperatures) or 
by using pulsed techniques^. The existence of phase 
slip centers/lines was confirmed in high-temperature 
superconductors^ and they were found to lead to S- 
shaped I-V characteristics in the voltage driven regime^. 



Recently phase slip lines were directly observed^ in low 
temperature superconducting stripes. In Refsi 12 ' 13 a new 
type of vortex dynamics (so called vortex channelling^ or 
the appearance of " kinematical" vortices^) was proposed 
which, in our opinion, is the "bridge" between slow vortex 
motion and the fast phase slip line regime. Furthermore, 
in Refi^ the experimental observation of quantum phase 
slips was claimed. 




FIG. 1: SEM micrograph of an Al square with lateral dimen- 
sion of 2 /Ltm. 

In this paper we present experimental results on the 
current-voltage characteristics of a mesoscopic supercon- 
ducting square with leads. In this geometry, the contacts 
play a crucial role (see Fig. 1) because the current den- 
sity is maximal in them. Therefore we cannot consider 
those contacts in equilibrium. This is essentially differ- 
ent from previous work on superconducting film/bridge 
attached to " massive" superconducting " banks" where it 
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was possible to assume the contacts in equilibrium. At 
low magnetic fields phase slip centers will appear in our 
sample in the narrowest places, where the current density 
is maximal. Because of inevitable heat dissipation the 
sample can transit locally to the normal state (for weak 
heat transfer) or to the superconducting resistive state 
(for strong heat transfer). In both cases the I-V charac- 
teristics are hysteretic due to heat dissipation and/or the 
dynamics of the superconducting condensate. Additional 
complications come from the effect of the magnetic field 
induced currents in the square on the phase slip process 
in the contacts. This makes our system new and to our 
knowledge this situation was not studied before either 
experimentally or theoretically. 

The paper is organized as follows. In Sec. II we present 
our experimental results and in Sec. Ill we give their 
interpretation on the basis of a solution of the time- 
dependent Ginzburg-Landau equations coupled with the 
equation for heat dissipation. Our conclusions are given 
in Sec. V. 



II. EXPERIMENT 

In Fig. ^ a SEM micrograph of our superconducting 
square made from Al using e-beam lithography is shown. 
The coherence length determined from a macroscopic co- 
evaporated sample was found to be £(0) = 156 nm. The 
thickness was 39 nm found from AFM and from X-ray 
measurements. Wedge shaped contacts with an opening 
angle of T = 15° were used. This shape was used to min- 
imize the effect of the contacts on the superconducting 
properties of the squar o 15 ; 16 . 

The experimental I — V characteristics are obtained 
by superimposing a small ac current (0.1 fiA rms) to a 
dc current Idc- The ac differential resistance is measured 
with a EG&G PAR 124A lock-in amplifier. The dc cur- 
rent is swept from negative to positive value. In order to 
ensure that the sample is in the normal state a high dc 
current is sent through the sample prior to the current 
sweep. Such sweeps are repeated for different magnetic 
fields. 

In Fig.[5]the differential resistance dV/dI(Id c ) is shown 
for four values of the magnetic fields. A clear hysteretic 
behavior is observed. When starting in the normal state 
and decreasing the current, the sample remains in a re- 
sistive state up to low currents. This resistive state is not 
the normal state anymore since a non-constant differen- 
tial resistance is observed. When starting from the super- 
conducting state and gradually increasing the current a 
non-resistive state is observed up to high currents. Con- 
trary to the transition seen at negative currents, a sharp 
transition from the non-resistive to the normal state is 
measured. While in the negative part the transition to 
the non-resistive state is accompanied by the appearance 
of a sharp peak, this is not observed for positive currents. 
This can be explained by our measuring technique and 
by the observed hysteretic behavior. When the transition 



occurs, the sample remains in the resistive state even 
when decreasing slightly the current so that the mea- 
sured ac voltage will either reflect the non-resistive or 
the resistive state, but not the transition. Above 1.2 mT, 
when the first vortex enters the sample, the transition to 
the resistive state is preceded by a small increase of the 
differential resistance [see Fig. |2£b)]. The shape of this 
part strongly depends on the vorticity of the sample sug- 
gesting a dissipation caused by vortex motion. At high 
magnetic fields, [see arrows in Fig. |2Jb)] small symmet- 
ric features are observed at high positive and negative 
currents. 

In order to study the magnetic field dependence in 
more details a contour plot of the differential resistance 
is shown in Fig. [3] as a function of the dc current and 
the magnetic field. The uniform gray (blue) areas at the 
left and the right correspond to the normal state and 
the white area in the middle to the non-resistive state. 
Oscillations are observed for the transition to the resis- 
tive state for positive currents with cusps each time the 
vorticity is changed [see arrows in Fig. b)] . The dissi- 
pation caused by vortex motion, which was already dis- 
cussed, can also be seen from the contour plot [light gray 
(green) area in Fig. [3] . The dissipation is only observed 
for a finite vorticity. At the transition from vorticity 
L=l to L=2 at /zo-ff=2.05 mT, the onset of the dissipa- 
tion seems to be continuous, but the dissipation increases 
more rapidly with increasing current when increasing the 
vorticity as can also be seen from Fig. Hfb). 

Beginning at ±2.5 mT small symmetric features are 
seen at high current [white (light blue) lines and the 
dark gray (dark blue) area in Figs. Eta)-(b)]. These fea- 
tures are visible up to ±20 mT and can thus not be aris- 
ing from the square since the square has a critical field 
H* 3 — 6 mT. It can only be generated by the contacts 
which have a higher critical field since a superconduct- 
ing/normal boundary with a sharp angle increases the 
critical fieldii±£. 

Above ±4.5 mT no hysteretic behavior is seen. All 
these features are perfectly reproduced on a different 
sample and are not dependent on the measuring condi- 
tions (sweep rate, amplitude of ac current, temperature). 
The same behavior is seen when fixing the applied dc cur- 
rent and sweeping the magnetic field: when the square 
is in the normal state, the sample remains in a resis- 
tive state up to low magnetic fields while when starting 
from the non-resistive state and increasing the magnetic 
field the resistive state is reached only at higher magnetic 
fields. 



III. THEORY 

To understand the experimental results we studied 
the current-voltage characteristics of two-dimensional 
superconductors using the generalized time-dependent 
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FIG. 2: Measured dV/dl as a function of the dc current at (a) noH=0 mT, (b) (m>H=3.2 mT, (c) fj,oH=4.2 mT and (d) (jloH=5 
mT measured at T = 1.100 K. The full black and dashed blue lines are for increasing and decreasing current, respectively. 



Ginzburg-Landau (TDGL) equation 19 

(d . 7 2 ^| 2 \ . 

7TWF U +ty+ T-arJ^ sa 

= (v-ia)V + (i-t-m 2 )^ (i) 

where the parameter 7 = 2r^;Ao/S is the product of the 
inelastic collision time te for electron-phonon scatter- 
ing and A = 4^?^ it 1 / 2 /n is the value of the gap at 
T = which follows from Gor'kov's derivation 2 ^ of the 
Ginzburg-Landau equations. 

This equation should be supplemented with the equa- 
tion for the electrostatic potential 

Aip = div(Im(V>*(V-iA)V>)), (2) 

which is nothing else than the condition for the conser- 
vation of the total current in the wire, i.e. divj = 0. In 
Eqs. (1,2) all the physical quantities (order parameter 
tp = \xj}\e 1 ^ , vector potential A and electrostatical po- 
tential tp) are measured in dimensionless units: temper- 
ature in units of the critical temperature T c , the vector 
potential A and the momentum of the superconducting 
condensate p = V0 — A are scaled in units $o/(27r£(0)) 
(where $0 is the quantum of magnetic flux) , the order pa- 
rameter in units of Ao and the coordinates are in units of 
the coherence length £(0) = (8fc B T c /7r7i.D) _1/2 . In these 
units the magnetic field is scaled with H C 2 = <I>o/27r£(0) 2 
and the current density with jo = cr„fi,/2erGL(0)£(0). 
Time is scaled in units of the Ginzburg-Landau relax- 
ation time rci(0) = Trh/8kBT c u, the electrostatic poten- 
tial (tp) is in units of ipo = h/2eTQL(0) (a n is the normal- 
state conductivity, and D is the diffusion constant). The 
parameter u is equal to 5.79 in accordance with Ref. 19 
and we used 7 = 40. We put A = (Hx,0,0) in Eqs. 
(1,2) because we limit ourselves to the case when the ef- 
fect of the self-induced magnetic field is negligible. This 
is valid in the experimental situation because the width 
of the sample is much less than the characteristic length 
A = A(0) 2 jdf (df is the thickness of the sample). 



Strictly speaking Eq. (1) is valid only very close to 
the critical temperature (see estimates for different low- 
temperature superconductors in Ref. 2 —). For example for 
bulk 'clean' Al the validity of Eq. (1) was derived only 
for the range AT ~ 10~ 4 K near T c . However our Al sam- 
ples are in the 'dirty' limit due to the small value of the 
mean path length I. As follows from RefsiS 2 * 2 ^ the rela- 
tion between current density j, the absolute value of the 
order parameter \ip\ and the momentum p are quite close 
to the Ginzburg-Landau relation j — p(l — p 2 ) = \ip\ 2 p 
even when T— > for such samples. Besides, when we 
turn on the magnetic field and/or the transport current 
the density of states of quasi-particles differs from the 
Bardeen-Schriffer-Cooper dependence 2 ^ and can become 
gapless 2 ^^ for high enough magnetic fields and/or trans- 
port currents. In this case Eq. (1) should be valid at any 
temperature T < T c because they were actually derived 
in the gapless limit (with 7 = O) 2 ^ or for small value of 
the gap A(T) < k B T^. 

The actual value of 7 for Al should be about 10 3 be- 
cause the time te ~ 10 _8 s is quite large in this material. 
However, the use of such a large 7-value is important if 
we intend to compare quantitatively the theoretical and 
experimental values for the critical current. As will be 
shown below there exist two critical currents which we 
call first I c i and second I C 2 critical currents. The first 
critical current is the current at which the sample goes to 
the non-resistive state and its value strongly depends on 
the value of the parameter 7 (in case of strong heat dissi- 
pation - see text below) . The second critical current has 
the meaning of the current at which the superconducting 
state becomes unstable and it could be determined from 
a stability analysis of the stationary Ginzburg-Landau 
equations. It implies that the current l c i does not de- 
pend on the 7 value which is the main reason why we are 
able to find quantitative agreement between theory and 
experiment for the position of the cusps in the l c i{iV) 
dependence (see below). 

In our theoretical model we considered the geometry 
depicted in Fig. 4 which simulates the real experimen- 
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FIG. 3: (Color online) (a) measured differential resistance 
dV/dl as a function of the dc current for different magnetic 
fields measured at T = 1.100 K. The curves for non-zero field 
are shifted for clarity. A color map surface is given for com- 
parison with (b). (b) contour plot of the differential resistance 
dV/dl as a function of dc current and magnetic field measured 
at T = 1.100 K for low magnetic fields. The arrows indicate 
the points where the vorticity changes. 



tal samples. To simplify our model we used linear con- 
tacts instead of the wedge shaped contacts used in the 
experiment. The main difference between them is that 
for wedge shaped contacts the order parameter is more 
suppressed in point A at low magnetic fields because the 
current density in this region is maximal (see Fig. 1). In 
order to inject the current in our system we used normal 
metal-superconductor boundary conditions at the end of 
the leads, i.e. ip = and — Vip — j. At the other bound- 
aries we used the usual insulator-superconductor bound- 
ary conditions: (iV — A)i(j\ n = and V<p\ n = 0. 

We also took into account the change of the local tem- 
perature in the sample in the resistive state by adding 
the heat diffusion equation to Eqs. (1,2) 



dT 

GeU ~dt = K eff AT + d fjn/Vn-h(T-T ), (3) 



where Tq is the bath temperature, C e // = (D s C s + dfCf) 
is the effective heat capacity, K e ff = (D s k s + dfkf) is 
the effective heat conductivity coefficient, and the heat 
transfer coefficient h = k s /D s governs the heat removal 
from the sample. Here we used a model for the tempera- 
ture distribution in thin superconducting films discussed 
in details in Refc2£ and C s , Cf, k s , kf are the heat capac- 
ity and the heat conductivity of the substrate (subscript 
s) and film/sample (subscript f), respectively. In this 
model it is assumed that the thickness of the substrate 
and the film D s + df is much smaller than the healing 
length A h = y/K eff /h > D s + d f . 

If heat removal is strong enough (large value for the 
coefficient h) we can neglect the effects due to the lo- 
cal change of the temperature. In the opposite case the 
results will quantitatively depend on the ratio between 
the healing length and the sample parameters (width of 
the square and the ratio between the value of the current 
density in the contacts and in the square). We chose our 
parameters in such a way that it optimizes the calcula- 
tion time (small value of C e ff) and we considered cases 
of large, intermediate and small value of the coefficient 
h. In dimcnsionless units [the same as Eq. (1)] Eq. (3) 
may be written as follows 



C eff 



dT 



K eff AT + j 2 n -h(T-T ), (4) 



where C eff = (D s C s /d f + C f )T c a n /T GL (0)ji, K eff = 
(D s k s /d f + k f )T c a n /£ 2 {0)jl h = hT c cr n /d f j 2 and the 
temperature is measured in units of T c . If D s C s /df «C/ 
and D s k s /df « we can use the Wiedemann- Franz 
law as an estimate for Cf and kf and we obtain for 

C7ff = tt 3 /48 ~ 0.65 and % = tt 4 /A8u 2 ~ 0.06 at 
a temperature close to T c . These values should be con- 
sidered only as a very rough estimate for the real mag- 
nitudes because normally the following inequalities are 
valid: D s C s /df 3> Cf and D s k s /df ^> kf. Because 
of the uncertainty in the actual values of C s and k s we 
used the following values: C e ff = 0.03 (to optimize cal- 
culation time), K e ff = 0.06, and h — 2 ■ 10~ 3 (which 
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corresponds almost to full heat removal at T = 0.9), h — 
2-10~ 4 (intermediate heat removal) and h — 2-10~ 5 (weak 
heat removal)^ and a bath temperature of To = 0.9. 
We checked that our results only weakly depend on our 
choice of C e t / and K e ff. As a boundary condition to Eq. 
(4) we take VT|„ = which means that heat is mainly 
transferred to the substrate. Only at the boundary be- 
tween the normal metal and the superconductor we used 
boundary conditions with fixed temperature Tjvs = T . 
The healing length is equal to A h ~ llf(T = 0.9T C ) for 
h = 2 • 10~ 5 and is comparable to the size of the sample. 
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we mean H c2 (T = 0.9T C ). 

First we studied the behavior of our model geometry 
sample in an applied magnetic field with zero transport 
current. In Fig. 5 we present the dependence of the free 
energy of the square and the value of the order parame- 
ter in the center of the edge (i.e. at point C in Fig. 4) 
as a function of the magnetic field. It turned out that 
for the chosen parameters (width of the square is 6£ and 
the size of the wires is 24£ x £ which are close to experi- 
mental values with £ ~ 0.333/im and H C 2 ~ 2.95 mT) no 
single quantum vortex state exists in the square and only 
surface superconductivity nucleates, i.e. the giant vortex 
state is present at H>0.8H C 2. Superconductivity vanishes 
in the superconducting square at H > 2.71-ff C 2 which is 
much larger than the third critical field H c j, = 1.69i/ C 2 of 
a flat infinite surface, which is a consequence of the shape 
and finite size of our system. By this value the vorticity 
in the system is equal to 10 and superconductivity may 
survive in the contacts up to much higher values of the 
magnetic field. 




FIG. 4: Model geometry. Between points A and B phase 
slip centers appear at H< H C 2. In point C we trace out the 
dependence of on the magnetic field (see text below and 
Fig. 5). 

Although in our numerical calculations we used units 
normalized at T=0 we will discuss here mainly the situ- 
ation for a bath temperature of Tq = 0.9T C . Therefore, 
it is more convenient to express the different quantities 
in units normalized at T = Tq. For example, under H C 2 



FIG. 5: Calculated dependence of the free energy (of the 
ground state) and the order parameter in the center of the 
edge of the superconducting square. 

What will occur when we switch the transport current 
on? Let us consider first the situation when the heat 
removal is quite good and the temperature of the sam- 
ple is equal to the bath temperature. Then up to some 
magnetic field H*(H* ~ 0.92/f C 2 for our parameters) the 
transition to the resistive state from the superconducting 
state occurs via the appearance of phase slip centers in 
the contacts [see Figs. 6(a,b)]. Current I c % at which this 
occurs slightly depends on the applied field because the 
order parameter in the square depends on H (see Fig. 
5) due to the induced screening currents. Via proximity 
effects the variation in the order parameter in the square 
influences the order parameter in the contacts and hence 
the critical current for nucleation of the phase slip cen- 
ters. 

At higher magnetic fields the vortex flow regime starts 
at I = I C 2(H). At fields close to H* this regime switches 
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FIG. 6: Snap-shots of the order parameter distribution (dark 
color corresponds to the maximal value of ^ and grey to the 
minimal one) in a superconducting square with attached leads 
at different magnetic fields and I > J C 2(H). 



FIG. 7: Current-voltage characteristics of the superconduct- 
ing square with contacts calculated for different heat removal 
coefficients h and magnetic fields. Current is normalized by 
the value jodw, where d and w are the lead thickness and 
width, respectively. 



to the phase slip process near the contacts [see Fig. 6(c)] 
while at fields higher than H* with further increase of 
the current it switches to the phase slip line regime in 
the square [or more exactly a line along which the vor- 
tices move very fas1ii£ - see Fig. 6(d)]. It occurs because 
the order parameter is strongly suppressed in the square 
by the magnetic field while in the contacts its influence is 
less pronounced due to the small width of the contacts. 
The voltage exhibits a jump at the transition from slow 
flux flow to the phase slip line/center regime. The larger 
H the smaller this jump in the voltagei 3 -. The actual 
value of H* , at which this change in the mechanism of 
destruction of superconducting state at / = I C 2(H) oc- 
curs, depends on the width of the contact. The narrower 
the contacts the higher the field at which flux flow in the 
square starts before the phase slip process occurs in the 
contacts. 

We should note here that the position of the phase 
slip line in the contacts depends on the applied magnetic 
field. When the order parameter at the edge of the square 
decreases, the phase slip center approaches the square 
and vice versa. For the case of wedge contacts we do 
not expect such a behavior because in this case the order 
parameter is more suppressed around the narrowest point 
where the current density is maximal. 

In Fig. 7 we present the I-V characteristics calculated 



for different values of the magnetic field. For low mag- 
netic fields and large heat transfer coefficient h the hys- 
teresis in the I-V characteristics is an internal property 
of the phase slip proces a 4 ' 7 i 13 . At the current I c \ (H) the 
sample goes to the superconducting non-resistive state 
from the phase slip regime at fields less than some crit- 
ical value H** or from the flux flow regime for higher 
fields. The actual value of the field H** depends on the 
value of the coefficient 7 and the width of the contact 
(see the above discussion for field H*). With increasing 
7, the minimal current at which the phase slip process is 
still possible, decreases and hence the field H** increases, 
because the current at which the flux flow starts in the 
sample does not depend on the relaxation times of the 
superconducting condensate. 

The magnetic field effects the hysteresis through the lo- 
cal suppression of the order parameteri This is the main 
origin of the decreasing and finally the disappearance of 
the hysteresis at the transition from the superconducting 
state or slow vortex flow regime to the phase slip regime 
at high magnetic fields. Another effect of the magnetic 
field is the slow increase of l c \ at low magnetic fields. The 
reason for this is the same as was found in Ref. 13 - the 
nonuniform current density distribution in the contacts 
due to the applied magnetic field. 
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FIG. 8: Calculated dependence of the temperature of the 
superconducting square at different currents and magnetic 
fields. Results for h — 2 ■ 10 -3 are not presented due to 
quite small changes (less than 1%) in the temperature. 



actual temperature of the sample is still larger than To 
(see Fig. 8). Consequently the sample is in the resistive 
state with a resistance less than the normal one. The 
range of currents for which such a process is possible de- 
pends on many parameters. For example, it increases 
with increasing magnetic field and decreasing heat trans- 
fer coefficient. 

At high magnetic fields the critical currents l c i and I c \ 
are quite small and even in the case of weak heat removal 
the TV characteristics resemble the ones with strong heat 
removal at low currents [Fig. 7(c)] due to weak heating 
[see Fig. 8(c)]. 
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When the heat removal is not effective, then the TV 
characteristics have a different shape at low and interme- 
diate magnetic fields. In Fig. 7 we present our results 
for two relatively small values of the heat transfer coef- 
ficient h. At low and intermediate magnetic fields the 
whole sample goes into the normal state (with T > T c or 
T < T c ) at the current / = I C 2(H) because of the large 
heat dissipation which is connected with the phase slip 
process. The value of that critical current is the same for 
any heat transfer coefficient h due to the absence of heat 
dissipation in the "pure" superconducting state. 

When we decrease the current, the temperature in the 
sample can become less than T c (see Fig. 8) while the 
sample will not go into the superconducting state be- 
cause at this temperature the current in the sample is 
too high in order that superconductivity can sustain it. 
Only when the temperature in the sample becomes less 
than some critical temperature T* < T c (which depends 
on the value of the heat transfer coefficient h) then su- 
perconductivity starts to nucleate in the square in places 
where the current density is minimal (near the corners 
of the square). For this current, magnetic field, bath 
temperature and in the absence of local heating, the flux 
flow or phase slip processes are impossible (see Fig. 7 
for high value of h) and the sample should go to the 
superconducting state. But due to heat dissipation the 



FIG. 9: Calculated dependence of the critical currents I C 2 
and 7 c i (for h = 2 • 10~ 3 and h = 2 ■ 1CT 5 ) on the applied 
magnetic field. I C 2 and I c i coincide for H >3.9 mT (1.32iS" C 2) 
at strong heat transfer (h = 2 • 10 -3 ) and for H >4.6 mT 
(1.56 H C 2) at weak heat transfer (h = 2 ■ 10 ). There exists 
a nonmonotonous behavior both in the I C 2(H) and I c i(H) 
dependencies which is connected with a change in the vorticity 
in the superconducting square. 

And finally in Fig. 9 we present the dependence of 
the currents I c i and I c \ on the applied magnetic field for 
strong and weak heat removal. There is a good quantita- 
tive agreement between the position of the cusps in the 
experimental and the theoretical dependencies of I C 2{H). 
We explain these cusps by abrupt changes in the vorticity 
and hence sharp changes in the order parameter distri- 
bution (see Fig. 5). The main difference between theory 
and experiment is in the amplitude of the variation of 
I C 2 with H (the theory predicts larger values) and in the 
value of this critical current. We believe that this dis- 
agreement originates mainly from the difference in the 
real shape of the attached contacts (Fig. 1) and the con- 
tacts used in our model (Fig. 4). In the experimental 
case there is a voltage lead which 'strengthen' the su- 
perconducting property near the narrowest point (as it 
was shown in the quasi-one dimensional limit in Rcf. 30 ) 
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and actually shifts the position of the phase slip center 
further from the square. Therefore, the effect of the vari- 
ation of the order parameter in the square (with applied 
field) should be less pronounced on the phase slip pro- 
cess in the contacts. Variations in I c2 become stronger 
at H > H C 2 when the resistive regime starts from the 
vortex flow regime in the square and hence the effect of 
a change in the vorticity is more 'visible'. 

We interpret the black right line in Fig. 3 as the line 
corresponding to the transition from the vortex flow to 
the phase slip line regime induced by increasing current 
at which a jump in the heat dissipation occurs. Because 
the heat removal is not effective, the square will become 
normal. Otherwise, it will be in the superconducting 
resistive state with a resistance close to the normal one. 
Here we would like to stress the following. In aluminium 
the decay length of the charge imbalance Aq or, in other 
words, the region where the normal current density is 
finite near the phase slip center— is quite larga^i Aq ~ 
50/iTO. The size of our sample is much less than Aq. So 
in this case it is quite difficult to distinguish between the 
normal and the superconducting resistive state because 
the differential resistance would be the same and equal 
to the normal one (see Ref. [30] in Ref. 7 ). 

With decreasing current (in absolute value) we cross 
the left black line (Fig. 3) and at low magnetic fields 
almost immediately enter the zero resistance state. Ac- 
tually, the current I c \ even increases a little. The same 
increase of I c \ can be reproduced theoretically if we as- 
sume that these transitions occur due to the decay of the 
phase slip process in the contacts at strong heat removal 
(see Fig. 9). 

A comparison with the experiment also shows that 
with increasing H we cross the left black line and do not 
enter the zero resistance state. Thus, we are probably in 
the vortex flow state but with the sample temperature 
larger than the bath temperature. The sample can be 
in the resistive state up to lower currents than it can be 
at the bath temperature (see Figs. 7,8 for small heat 
transfer coefficient). 

At very high fields dissipation is not very important at 
currents close to the critical ones I = I C 2(H) = I c i{H) 
(due to their small value) and besides there is no 'inter- 
nal' hysteresis due to the phase slip process. As a result 
there is no hysteresis in the current when the non-zero 
resistance state appears and when a fast change in the 
resistance occurs in our sample. 

IV. CONCLUSION 

Hysteresis in the current-voltage characteristics of su- 
perconducting wires, films or mesoscopic samples may 
appear due to heat dissipation or/and due to 'inter- 
nal' hysteresis connected with the existence of phase slip 
lines/centers. We believe that in our measurements we 
have both types of hysteresis which are responsible for 



the observed effects. At low magnetic fields the sam- 
ple enters the normal state due to the appearance of the 
phase slip process and a strong heat dissipation at the 
critical current I C 2- Because of a discrete change of the 
vorticity in the superconducting square the order param- 
eter changes abruptly at some values of H and it leads 
to cusps in the dependence of I C 2(H). At higher fields, 
instead of a high dissipative phase slip process, we have 
slow vortex motion and heat dissipation results in a weak 
effect on the I-V characteristics for / > I c2 (H). 

When we decrease the current (at fixed value of the 
magnetic field) the sample goes first from the normal 
to the superconducting resistive state (left black line in 
Fig. 3) and than slow vortex flow starts in the supercon- 
ducting strip (at high magnetic fields) or the phase slip 
process in the contacts (at low magnetic field). Because 
the temperature of the sample may be higher than the 
bath temperature, the current, at which the sample goes 
to the nonresistive state, may be smaller than I C 2 even if 
the resistive state starts as a vortex flow. So actually we 
'need' heat dissipation to explain this effect. 

One of the main results of our paper is that we 
show (by a self-consistent solution of the time-dependent 
Ginzburg-Landau equation and the heat diffusion equa- 
tion) that heat dissipation does not necessarily lead to the 
destruction of superconductivity as it was supposed in a 
recent papers (in order to explain recent experiments 
on MoGe films (see references therein)). In these sam- 
ples the value of the coefficient h is quite small due to 
the small value of the normal conductivity and both crit- 
ical currents I c ± and I C 2 are rather large because of the 
absence of the external magnetic field. In general, the 
situation may be more complicated when heating does 
not destroy superconductivity. In such a case, heating 
leads to additional complexity in the dynamics of the su- 
perconducting condensate due to the local heating of the 
sample. 

Another result is, that by comparing the experimental 
and theoretical I C 2 (H) dependence we may distinguish 33 
between the giant vortex and the single quantum vortex 
state which appear in the sample. It allows, in princi- 
ple, to study experimentally the transformation between 
these two different configurations as a function of the 
shape, size of the sample and external magnetic field. 
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